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abstract 

We study how the initial inhomogeneities of the universe affect the 
onset of inflation in the closed universe. We consider the model 
of a chaotic inflation which is driven by a massive scalar field. In 
order to construct an inhomogeneous universe model, we use the 
long wavelength approximation ( the gradient expansion method ). 
We show the condition of the inhomogeneities for the universe to 
enter the inflationary phase. 
PACS number(s): 98.80.Cq, 98.80.Hw 



I. INTRODUCTION 

The main objective of inflationary scenario [Q is to explain why the present universe is homogeneous and isotropic 
on large scales without a fine tuning of initial conditions. It is important to study the generality of initial conditions 
for the onset of inflation. 

Many works on inflation assume homogeneous and flat universe at the pre-inflationary phase 0]. In this case, 
inflation is generic. The spatial curvature effect on the occurrence of inflation is investigated in the Friedmann- 
Robertoson- Walker model with a massive scalar field by Belinsky et al and in the Bianchi IX model with a 
cosmological constant by Wald Q. They concluded that inflation is not a general property and only the positive 
spatial curvature can prevent the universe from inflating. 

On the other hand, the role of initial inhomogeneities on the occurrence of inflation was studied by a linear 
perturbation analysis Q and by a numerical simulation with special symmetries . In ref . || , Goldwirth and Piran 
concluded that the crucial feature necessary for inflation is a sufficient high average value of the scalar field which 
drives inflation over a region of several Hubble radius. 

We study how the initial inhomogeneities of the universe affect the onset of inflation by use of an alternative 
approach. In order to treat the inhomogeneities of the universe, we use the long wavelength approximation. 

In the long wavelength approximation we take a synchronous reference frame where the line element is of the form 

ds 2 = g^dx^dx" = -dt 2 + ^(t, x)dx i dx j . (1.1) 

Throughout this paper Latin letters will denote spatial indices and Greek letters spacetime indices. In the usual 
method, we neglect all spatial gradients of jij (t, x) in the Einstein equations and construct an approximate solution 
with the characteristic scale of inhomogeneities larger than the Hubble radius Furthermore, we consider the 

quasi-isotropic universe as a special case in the form 

ds 2 = -dt 2 + a 2 {t)h lJ {x)dx l dx 1 , (1.2) 
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where hij(x) is an arbitrary function called 'seed metric'. The universe is assumed to be filled with a perfect fluid 
characterized by energy density p and pressure p = (T — I) p. In this case the Einstein equations reduce to 
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where a dot denotes the derivative with respect to t and n = 8ttG. These are nothing but the equations which the 
scale fa ctor a(t) of the flat Fried mann univ erse should satisfy. We call the approximate inhomogeneous solution given 
by Eq.( |1.2| ) satisfying Eqs.( |1.3D and (|1.4|) the locally flat Friedmann solution and it is used as a starting point to 
solve the Einstein equations iteratively |10| , pT| . In the iteration scheme, the spatial metric can be expanded as a sum 
of spatial tensors of increasing order in spatial gradients of hij(x). Thus the approximation is called the gradient 
expansion method. 

Recently, the influence of initial inhomogeneities on the occurrence of inflation is studied by using the gradient 
expansion method ]l3|-|l5|]. From investigations in the homogeneous universe model, as noted above, the effect of 
positive curvature is important because it can prevent inflation. So, we should investigate the inhomogeneities of 
universe with a non-small positive spatial curvature to clarify the generality of inflation. However, the lowest order 
solution in the gradient expansion method is a locally flat Friedmann universe and the spatial curvature is treated as 
a small quantity in the expansion scheme. 

In ref. [jl5|] , the gradient expansion method from the locally closed Friedmann universe is introduced and the 
inhomogeneous closed universe with a cosmological constant is studied as an application of the scheme. This method 
is useful to treat the inhomogeneities of non-small curvature associated with a three-space which is conformal to the 
constant curvature space at lowest order. In this paper we consider inflation which is driven by a massive scalar 
field in the inhomogeneous universe by use of the first order approximation of the gradient expansion from the locally 
closed Friedmann universe. We will show the condition of the inhomogeneities of the spatial curvature for the universe 
to enter the inflationary phase. 

This paper is organized as follows. In Sec. 2, we derive basic equations by use of the first order approximation in the 
gradient expansion method from the locally closed Friedmann universe. In Sec. 3, we show numerical results. Sec. 4 is 
devoted to conclusion. We use units in which c — Ti — 1. 



II. BASIC EQUATIONS IN IMPROVED GRADIENT EXPANSION METHOD 



We take the matter in the universe to be a scalar field minimally coupled to gravity. The energy-momentum tensor 
of a scalar field is given by 
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where V((f>) is the potential of the scalar field. In the synchronous reference frame Eq.( |l.l[) , the Einstein equations 
for the gravitational field coupled to the scalar field read 
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(2.2) 
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where 3 i?^ is the Ricci tensor associated with 7^, = Ty, K = Y J ^iji 7 = detTy, and a semicolon denotes the 
covariant derivative with respect to 7y . The equation of motion for the scalar field is 



dV(<P) 



= 0. 



(2.5) 



Now, we assume the metric in the form |15| 
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a 2 (t, £l(x))hij(x), 
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where hij(x) is the metric of thre e-di mensional sphere whose Ricci curvature is Rij(hki) = 2hij. The scale factor 



a(t, 9l(x)) in the spatial metric Eg. ( |2.q ) has the space-dependence by an arbitrary function Q(x). In the usual gradient 

(0) 

expansion method, all spatial gradients of 7 y (t, x) are assumed to be much smaller than time derivatives. Roughly 
speaking, the spatial curvature 3 R is assumed to be smaller than K 2 and it is neglected in the lowest approxima tion . 
Here, we consider slightly different approximation. The spatial scalar curvature for the metric in the form of Eq.( [2.6| ) 
is 
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where a prime denotes the derivative with respect to Q and V; denotes the covariant derivative with respect to h{j. 
We assume that the second term in the right hand side of Eq.(2.7) which contains spatial gradients of is smaller 



,(0),(0) 



than the first term. If we set ij(7y) = 6a (t, fl) Eq.(2.7) is rewritten in the form 



3 i?(7 lJ ) = 3 i?(7 y ) 
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Thus the case, we consider here, is that the spatial variation of the spatial curvature is smaller than the value of 
itself. All spatial gradi ents of il(x) is neglected in the first step of our approach. It should be noted that even if 
the second term in Eq.( |2.7| ) or fl2.8| ) exceeds the first term, the approximation reduces to the usual long wavelength 
approximation so far as these are smaller than K 2 . 

Though the metric in the form of Eq.(2.6) has not adequate degree of gravitational freedom, the simple form on 



which we concentrate is useful to investigate the inhomogeneities of the spatial curvature. It might be very complicate 
to treat the general form of inhomogeneities of curvature Jl2]| . 



F rom now on, we consider that the spatial metric in the lowest order are described by Eq.(2.6). Substituting 
Eq.( [2.6|) into Eq^^), and neglecting all spatial derivatives of Q(x), we obtain the equation which a and tfi must 
satisfy: 



d(t,n(x)) 
a(t,rt(x)) 



-4> 2 {t, x) + V(<l>(t,x)) 



(2.9) 



From Eq.(2.9), it is natural that the scalar field has the space-dependence through Q(x). We assume that the form 
of the scalar field in the lowest order is 

C ?(t,a:) = ^o(* i n(x)). (2.10) 
Substituting Eqs.(2.6) and ( p. 10 ) into Eqs.(2.3) and (2.5), and neglecting all spatial derivatives of fi(x), we have 

a(t,Q) fa(t,fl)\ 2 2 
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(2.11) 
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Equation ^. 4| ) is trivial in the lowest order. Equations (2^), ( 2.11 ) and ( 2.12 ) have the same form of the equations for 
the closed Friedmann universe with a homogeneous scalar field. 

At next order, we consider the second order in spatial gradients of Sl(x). We will take corrections to the metric and 
the scalar field of the form 



7y (t,x) = a 2 (t,il) 



9. 
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where ViVjil/fi and WiilWjft/fl 2 are traceless tensors denned by 
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Substituting 7y = 7 jj + 7 y and (f) =</)-{- </) given by E ns. (B. g ), ( |2.13| ), Q2.10| ) and fl2.14| ) into the Einstein equations 



and the equation of motion for the scalar field Eqs.(2.2)-(f2jj J7"& n d comparing the coefficients of the derivative of fi, 
we get the first order equations which govern the evolution of variables F, F, G, G, P and Q. We should perform 

The expression of it in the first order approximation 



tedious calculation to get 3 Rj associated with 7y 
appears in Appe ndix [A| 
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where ^— is the lowest order value of H = a/a, G = G + F 2 and G = G — \F 2 . Note readers should pay 
d(f> o dq> 

attention to the commutation relation of the covariant derivatives which is not vanishing because the spatial surface 
in the lowest order solution has a positive curvature. From Eq.(|2.5[), we have 
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From Eqs.(2.2) and (2.4) we obtain 



F = 2 



a a 
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r = - -<>nF + ^(j>' n, (2.i9) 

(2.20) 

n + K^ (j)' n, (2.2i) 

dV , 
F + 2HF = 2k -2<j> n P+ — P , (2.22) 

(2.23) 

There are nine equations (|2.15 )— ( 2.23D for the six variables F, F,G,G, P and Q. All of these equations arc not 

inde pende nt. Six equations govern the evolution and three are constraint conditions. So we find a solution of 

Eqs.(2.15)-(2.20) with an initial condition which satisfies the constraints. In order to solve these equations, we need 
to know the time evolution of a', a", <j>' n , <j>'^ and F . Differentiating Eqs.( |2.9D , (2.11), (2.12) and (2.21), we obtain the 
equations for them, Eqs.(B2)-( |B6| ), which appear in Appendix [b|. At higher order, we consider the corrections which 
are constructed by the spatial derivatives of Q(x) and solve the Einstein equations order by order in the gradient 
expansion. The spatial metric and the scalar field can be expanded as 
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where the notation (V( 2p )f2)^ denotes symbolically symmetric spatial tensors which contain 2p spatial gradients of 
the suffix A distinguishes the tensors belonging to the same class, and a(t, Q), (j>o(t, ft), FA p ){^ O) and P^it, fi) 
are function of t and Q(x) which is determined by the Einstein equations. In the case of p = 1, 2, the explicit form of 
(V (2p) fi)^ appears in the Appendix 

This expansion scheme is valid when the characteristic scale of inhomogeneities of the spatial curvature is larger 

than the Hubble radius (a/a) -1 or the spatial curvature radius 3 R ~ x l 2 . If the Hubble radius is smaller than the 
spatial curvature radius, this expansion scheme reduces to the usual gradient expansion. 



III. THE INFLUENCE OF INHOMOGENEITIES FOR INFLATION 



We consider a chaotic inflation model with a massive scalar field, i.e., V(<j>) = ^m 2 <p 2 as an inflaton field. 
When the energy of the scalar field is domi nated by a kinetic term at the early stage, the initial behavior of solution 
without inhomogeneities for Eqs.(|]^), ( 2.11 ) and ( 2.12| ) is expressed by series expansions in the form 
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(3.1) 
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where (f>o is the constant which corresponds to the freedom of initial value of the scalar field. The time variable t is 
non-dimensionalized by use of m. On assumption that we neglect all spatial derivatives of f2, we can generalize the 
solution Eqs.(3.1) and (3.2) in the form 
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(3.3) 
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The local scale factor and the scalar field are expressed by series expansions of (m~ 2 /Vt 2 )t^ in the initial stage. At the 
beginning of the universe, the local scale factor behaves as a(t, Q) oc t 1 ^ 3 and after that the time evolution depends 
on the local value of Q(x). In the initial stage where the kinetic energy of scalar field dominates the potential energy, 
the behavior of a(t, H) does not depend on 4>q but in the late stage where the potential energy becomes effective, it 
depends on </>o- 

We can see the time evolutions of these variables, which is characterized by the value of Q(x) and (fio, by solving 
Eqs.(2.9), ( 2.11 ) and ( 2.12 ) numerically. Fig.|l| shows the inflationary region and the recollapsing one on (il 2 , </>q) plane. 
The local scale factor with parameters in the inflationary region enters the accelerating expansion phase: a/a > and 
a I a > 0. On the other hand, the local scale factor with parameters in the recollapsing region enters the recollapsing 
phase: a/ a < and a /a < 0. From FigJU we see that inflation favors large f2 and large \<fro\ as is expected by the 
investigation of the homogeneous model pf. For fixed 0o there is a critical value of f2, f2 cr , if f2 > fl cr inflation occurs. 
Since the critical value il cr depends on <po then the condition for the onset of inflation in the lowest order solution is 
written in the form 



1 1 







(3.5) 



In the Fig.^ f2 cr (0o) is drawn as the boundaries between the inflationary region and the recollapsing region. 

In the initial stage where the local scale factor a(t, Cl(x)) is proportional to Q(x)t'S , the spatial scalar curvature 
behaves as £1~ 2 £~3 . The condition (pjj) is, then, translated to the condition on the initial curvature as 



3 Rinit(t,Q) < 3 R C r(t,(j>o), 



(3.6) 



where 3 R cr = 6t~ifl~ 2 ((j>o). This condition for the onset of inflation is a restriction on on l y the local value of 3 R. 

Next, we take the first order correction terms into account. We integrate Eqs. ( 2. 15 )-( |2.20 ) and Eqs.(B2)-(B6 ) 
numerically. We restrict ourselves to consider growing solution whose initial asymptotic behavior is 
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The time evolution of F, G, P and Q is shown in Figs.2-(i)-2-(iv). 

In the case that the local scale factor in the lowest order will enter the inflationary phase (ft 2 > ft^), F and G grow 
in proportion to time, and P and Q approach constant values. On the other hand, in the case that the local scale 
factor in the lowest order will enter the recollapsing phase (ft 2 < ftj! r ), F, G, P and Q diverge and the approximation 
breaks down in the course of time. These behavior are consistent with the re sults by the linear perturbation analysis 
||. The trajectory of scalar field in (0, </>) space is shown in Figs j3-(i)| - |3-(ii) . 

At each spatial point, we define an effective local scale factor a e jj(t,x) by 
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which describes how the universe expands at the spatial point. We can follow the evolution of a e jj(t, x) by the evolution 

of a(t, ft), F(t, ft) and G(t, ft). The behavior of a e Mt,x) is characterized by four parameters </>o,ft> V;V z ft/ft and 

VjftV'ft/ft 2 . In addition to the value of <fro and ft, the spatial derivatives of ft affect the evolution of the local scale 
factor. 

The effective local expansion rate and acceleration rate are given by 
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G(t,n) 



V/ftV'ft 1 
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We can divide the four-dimensional parameter space (</>o, ft, V/V'ft/ft, V/ftV'ft/ft 2 ) into two regions: inflationary 
region and recollapsing region by a three dimensional hyper surface. The effective local scale factor with the parameters 
in the inflationary region enters the accelerating expansion phase: d e jj/a e jj > and ^eff/ a eff ^ ^- O n the other 
hand, the effective local scale factor with the parameters in the recollapsing region enters the recollapsing phase: 
a e g/ a e jj < and a e jj/a e g< 0. We assume the approximation is valid while 
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(3.8) 



We evolve the universe numerically while both of these conditions Eq.(3.S) hold 



Evolving of the effective local expansion rate and acceleration rate numerically, we see how the two spatial gradients 
V;V ; ft/ft and V;ftV'ft /ft 2 , affect t he occurrence of inflation. For fixed </>o, the inflationary region and recollapsing 
region are shown in Figs 



OH-(iii) 



Between these two regions, there is a fuzzy region where error becomes large and 
one of inequality in Eq. ( |3.8| ) does not hold before the effective local scale factor enters the accelerating expansion phase 
or the recollapsing phase. We see that the positive (negative) VjV'ft/ft helps the universe to enter the inflationary 
(recollapsing) phase. On the other hand, V;ftV'ft/ft 2 tends to prevent the onset of inflation. 

From the numerical results, we obtain the time independent condition for the onset of inflation in the vicinity of 
ft/ft cr — 1 = and V;V'ft/ft = V;ftV'ft/ft 2 = in the parameter space, within the first order of gradient expansion, 
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(3.9) 



where a and j3 are constants. The parameter a shows the strength of the effect of inhomogeneities on the occurrence 
of inflation. In the massive scalar inflaton model, the value of a is about 2.0. The value of (3 depends on the initial 
value of the scalar field (/>o: P = 4.3 for <pa — 5.0k -1 / 2 and f3 = 3.3 for <fio — 10. Ok -1 / 2 . In the early stage where 
a 2 (t, ft) = ft 2 £s , the spatial curvature 3 Ri n it in the first order approximation is 
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By use of Eq.( |3.10| ), the criterion ( |3.9| ) is rewritten as 
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where 3 R cr (t,4>o) is the marginal value of the spatial curvature for inflation in the closed Friedmann model appears 
in Eq.( |3.6| ) and the parameters a and /3 are 

12a - 3ap - 7 



= (3a-2), $ = 



2(3a-2) 



(3.12) 



The criterion ( 3.11 ) is a condition on the initial value of the spatial curvature and its spatial derivatives. Because a 
and P are several numbers when <fto = (5.0 ~ 10. 0)k -1 / 2 , the onset of inflation is occur if the initial value of spatial 
curvature 3 Ri n it is less than the critical value 3 R cr (t,4>o) over a region which has several size of the local curvature 
radius ( 3 Rinit)~ 1/2 • 

Goldwirth and Piran studied how inhomogeneities influence the inflationary epoch numerically in the case of 
spherical symmetry ||. They conclude that the crucial feature necessary for inflation is a sufficiently high average 
value of the scalar field ( suitable value for inflation in homogeneous universe ) over a region of several Hubble radius. 
In order to compare their results with ours, we impose fl 2 (x) to spherically symmetric form, fi 2 (x ; 0, ip) = fi 2 (x), in 
the spherical coordinate where hij = diag[l, sin 2 sin 2 %sin 2 6*]. We consider fi 2 (x) in the form 



n 2 (x) = nl + (n 2 -nl) 
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exp [— A" 



2 cos 2 I] 
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where fig and fi 2 are the value of fi at \ — and \ — n i respectively, and A is a parameter. We fix 



7T, respectively, ana A is a parameter. We tix (po = 5.0k 
for example, then fl 2 cr « 0.70m~ 2 (see Fig§). The form of tt 2 for (fig + fi 2 )/2 = fl 2 cr and fig - fi 2 = 4 x lO^rrr 
shown in Fig.||. Varying A, we show inflationary and recollapsing regions in the three-dimensional space (see Fig.[ 
At the origin \ = universe can enter the inflationary phase whenever A is smaller than about 0.95. From Figl 
we see that the region of a suitable value of fi 2 is < \ < 1.85. Comparing the physical length L ~ a(t, fi)x for this 
region and the initial curvature radius (^Rmit)^ 1 ~ o(t, f2)/\/6, we get L ~ 4.5( 3 Ri„it) -1 / 2 . In this case, inflation 
occurs at the origin when the initial value of spatial curvature is less than the critical value over the region 4.5 times 
initial curvature radius. 



IV. CONCLUSION 



We consider a chaotic inflation model which is driven by a massive scalar field. In this model, we study how the 
initial inhomogeneities of the universe affect the onset of inflation by use of the gradient expansion from a locally 
closed Friedmann universe. 

At lowest order of our approximation, inhomogeneities of the metric and the scalar field are described by an 
arbitrary function fi(a;). At next order, using the effective local scale factor a £ jj{t,x), we investigate the effect of 
inhomogeneities of the universe for the onset of inflation. 

From numerical results, we obtain the condition for the onset of inflation. The spatial curvature in the early stage 
should be less than the critical curvature value 3 R cr (t, <j>o) over a region which has several times size of the local 
curvature radius ( 3 i?i ra it(i)) _1 ^ 2 for the onset of inflation in the future. 

We can compare the effect of inhomogeneities in the inflation model which is driven by a massive scalar field with 
the one in the model driven by a cosmological constant . The value of a in the scalar field model is larger than 
the one in the cosmological constant model. It means that the inflationary model driven by a massive scalar field is 
more sensitive to the inhomogeneities than one by a cosmological constant. 
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APPENDIX A: EXPRESSION FOR THE SPATIAL CURVATURE IN THE FIRST ORDER 



In the gradient expansion from the locally closed Friedmann universe, the spatial metric up to the second order of 
the spatial gradients is given in the form 
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The spatial Ricci curvature for the metric is 
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APPENDIX B: EQUATIONS FOR a , a" 



We derive the equati on fo r a' , a", </>q, <Pq and F 
From Eqs.(2.9) and (2.11) we obtain 
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Differentiating Eqs.( [2.12] ) and (Bl) with respect to f2, we find the equations which a'(t, Q.) and <p'(t, SI) must satisfy 
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Differentiating Eqs.(]B2|) and ( E3 ) again, we find 
.a . ,, I'd 2 + 1 k • 



a 

-2ko 

_ a 



a' 2 a'da' a' 2 (a 2 + 1) 
- + 2 — i— = - 



a " + 2 — - 4 



dV 



■, aa - aa ; 

D0o 5 O0 O 



a'd' a"d „da'' 



Similarly from Eq.(2.21), we find that F must satisfy 



— ' aa — a a „ 
F = 2 - + 2tt 



a" 
2- 



(B2) 
(B3) 



(B4) 
(B5) 



(B6) 
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APPENDIX C: THE SPATIAL TENSORS 



The spatial tensors (V^ 2p ^f2)^- which contain 2p spatial gradients of Q are given as follows. 
For (V< 2 )0)4 

__ and ____ 

For (V< 4 )0)4 

ViVj-nViV'n VjfiVjViV'n V;V»VjOv'n v.v^v^v^) 
VifiVjQViV'n ViV 3 -nvjnv J n Vifiv^V/fiv'n) 
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Figure Captions 

Figf _ 

In the lowest order approximation, an inflationary region and recollapsing region are shown in (Q 2 ,(j)o) plane. 



Fi g s. ra - p^) 



Typical examples of the time evolution of (i)F(t, 0), (ii)G(t, f2), (iii)P(t, f2) and (iv)Q(i, £1). 



<) 



The trajectories in (</>, </>) space. 

(i) ViOV'n/fi 2 = and VjV'fi/fi = 0.01(a), 0.001(b), 0(c), -0.001(d) and -0.01(e). 

(ii) V;V'rj/0 = and V^V^/fi 2 = 0(a), 0.0001(b), 0.001(c) and 0.01(d). 
At the endpoint of the trajectory the approximation breaks down. 



Figs. F(i)l - Riii)| 

An inflationary region and a recollapsing region for fixed 0o = 5.0k- 1 / 2 {fl 2 cr w 0.70m- 2 ) are shown in the parameter 
space ( V;V'f2/J7, V;57V'J7/f2 2 ). A meshed region is a region where the approximation breaks down. 



Fig.|. 

The shape of 51 2 (x) f° r varying A. 
Figf 

An inflationary region and a recollapsing region in the three-dimensional space for fixed </>o 
the case that fl 2 is given in Fig.|[ 



5.0k 1 I 2 are shown in 
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n 2 [m ! ] 

Fig. 1. 
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Fig. 2-(i). 
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Fig. 2-(ii). 
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Fig. 2-(iii). 
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Fig. 2-(iv). 
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Fig. 3-(i). 
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Fig. 3-(ii). 
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Fig. 4-(iii). 
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Fig. 5. 
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Fig. 6. 
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